Finite temperature spin-dynamics and phase transitions in spin-orbital models 



0\ ■ 
O ■ 

o ■ 
(N : 

> : 
o . 

<N : 



C.-C. Chen/'2 B. Moritz,- J. van den Brink,--^ T. P. Devereaux,' ^ and R. R. R Singh"* * 

' Department of Physics and Geballe Laboratory for Advanced Materials, Stanford University, CA 94305, USA 
^Stanford Institute for Materials and Energy Science, SLAC National 
Accelerator Laboratory, 2575 Sand Hill Road, Menlo Park, CA 94025, USA 
^Institute Lorentzfor Theoretical Physics, Leiden University, P.O. Box 9506, 2300 RA Leiden, The Netherlands 
''Department of Physics, University of California, Davis, CA 95616, USA 
(Dated: November 2, 2009) 

We study finite temperature properties of a generic spin-orbital model relevant to transition metal compounds, 
having coupled quantum Heisenberg-spin and Ising-orbital degrees of freedom. The model system undergoes a 
phase transition, consistent with that of a 2D Ising model, to an orbitally ordered state at a temperature set by 
short-range magnetic order. At low temperatures the orbital degrees of freedom freeze-out and the model maps 
onto a quantum Heisenberg model. The onset of orbital excitations causes a rapid scrambling of the spin spectral 
weight away from coherent spin-waves, which leads to a sharp increase in uniform magnetic susceptibility just 
below the phase transition, reminiscent of the observed behavior in the Fe-pnictide materials. 

PACS numbers: 74.70.Dd, 75.10.Jm, 75.40.Cx, 75.40.Gb 
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Correlated materials exhibit intriguing phenomena arising 
from the interplay between spin, charge, lattice, and orbital 
degrees of freedom. Orbital degrees of freedom can emerge in 
multi-band systems such as 3d transition metal compounds. In 
these systems, spins and orbitals are strongly coupled as spin 
exchange is the dominant interaction between different orbital 
occupations, which in turn support different spin order This 
correlation can lead to a phase transition in one or both vari- 
ables, the collective effects of which can be antecedent or sub- 
sequent to a lattice structural transition [1]. A paradigmatic 
example is manganites where orbital ordering is essential in 
explaining the magnetic properties and phase transitions [2]. 

The newly discovered Fe-pnictide superconductors [3] dis- 
play superconductivity in close proximity to magnetic order 
The observed collinear (tt, 0) magnetic order [4] has been 
studied theoretically from both weak- [5] and strong-coupling 
points of view [6]. In particular, an anti-ferromagnetic (AF) 
coupled J1-J2 Heisenberg model on a 2D square lattice (de- 
picted in Fig. 1 (a)) can give rise to an AF (tt, 0) order when 
J2 > [7]. Alternatively, this (tt, 0) order may be ob- 

tained through an anisotropic Jia-Jib-J2 model [8], where 
one has strong AF coupling in the x-direction, and ferro- 
magnetic coupling along the y-direction, as shown Fig. 1 
(b). Interestingly, recent neutron scattering data [9] indicate 
that the magnon energy is a maximum at momentum transfer 
(tt, tt). This strongly favors the Jia-Jib-J2 scenario which 
reproduces the observed spin wave dispersion, see Fig. 1 (d). 

A possible microscopic origin for the anisotropy in the Jia- 
Jib-J2 model is orbital ordering [10, 11]. When the orbitals 
are ordered, the lattice distorts and the orbital lobe orientations 
can cause a vanishing effective hopping in certain directions, 
as in ID edge-sharing copper oxides [12]. In conjunction with 
double exchange [13], orbital ordering, affecting the direct 
and super-exchange processes, can lead to even sign-changing 
anisotropic exchange interactions. Proposals have been put 
forth that consider ordering between the Fe 3dxz and 3dyz or- 
bitals as a possible mechanism for the observed magnetism of 
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FIG. 1 : (Color online) Schematic representation of various Hamil- 
tonians: (a) The J1-J2, (b) the Jia-Jib-J2, and (c) the spin-orbital 
models. At zero temperature the additional orbital degrees of free- 
dom freeze-out, and the spin-orbital model reduces to the J\a-Jib-J2 
model, (d) The magnon dispersion cjr (in units of Ji / J\a) calculated 
from linear spin wave theory. Here Jib — — O.lJia, J2 ~ 0.4Jia in 
the J la- J lb- J 2 model, and J 2 ~ Ji in the J1-J2 model. The spin 
wave energy forms a maximum at (tt, tt) in the Jia- Jit- J2 model, 
which is a minimum in the J1-J2 model. 



the Fe pnictides [10, 11, 14]. These proposals remain contro- 
versial in part because early band structure calculations [8], 
which agree well with a variety of experiments, show a very 
small difference in the occupation of dxz and dyz orbitals in 
the magnetically ordered tetragonal calculation [15]. On the 
other hand, recent ab initio calculations suggest robust orbital 
order using Wannier orbitals [16]. Indeed, if the magnon en- 
ergy is a maximum at (tt, tt) as reported in Ref. [9], this im- 
plies not just a small anisotropy due to for example structural 
considerations, but an extreme sign-changing one associated 
with additional broken symmetry [11]. 

In this paper, we address the question: If the anisotropy in 
exchange constants observed in neutron scattering is related to 



2 



orbital order, what other consequences follow. To answer this, 
we consider the following spin-orbital Hamiltonian relevant to 
the schematic in Fig. 1 (c): 



H = Ji ^[Si ■ Sj+£ niUi+i + Si ■ Si+y (1 - ni){l - rii+y)] C 

i 



=-0.5 



<ij> 



(1) 



where is an Ising variable taking values or 1, and Si is a 
spin-i operator: ■ Si = S{S + 1), and S = ^. The sums 
<> and <C^ run over nearest- and second nearest-neighbors, 
respectively. This model describes a system consisting of two 
orbitals per site, with the occupation controlled by the Ising 
variables: ??i = represents orbital 1 occupied; rii = 1 repre- 
sents orbital 2 occupied. 

When the interactions are dominated by an AF coupled Ji , 
the above model finds its lowest energy configuration in a per- 
fect ferro-orbitally ordered state corresponding to all ni = 
or 1 . Therefore, at zero temperature T = this Hamiltonian 
reduces to the Jia-Jib-J2 model with Jia = Ji — J/ and 
Jit, ^ —Jf. On the other hand, the finite temperature proper- 
ties would be quite different due to orbital fluctuations and ex- 
citations. In the following calculations we take the parameters 
from neutron scattering data on CaFe2As2 [9]: SJia = 50 
meV, SJf = 6 meV, and SJ2 = 20 meV. 

We are interested in the finite temperature spin dynamics 
of these systems; however, there are few numerical methods 
capable of accomplishing this in a controlled manner. We use 
the exact diagonalization (ED) technique, which has been uti- 
lized extensively to investigate both zero and finite tempera- 
ture properties for various quantum lattice models [17]. We 
use N=16 site square plaquettes, already requiring a large 
computational effort due to the additional orbital degrees of 
freedom. Lattice translation, rotation, reflection and Ising- 
orbital inversion symmetries reduce the 2^^ Ising configura- 
tions to 733 distinct ones. We fully diagonalize the Hamilto- 
nian in these Ising sectors and calculate dynamical quantities. 

Our main results are as follows: (i) In a purely 2D system, 
where in accord with the Mermin- Wagner theorem the spin- 
rotational symmetry can not be spontaneously broken except 
at T = 0, the orbital degrees of freedom undergo a phase 
transition at a temperature scale ~ 0.2 Ji set by short range 
magnetic order, (ii) At temperatures below O.lJi, the Ising 
variables are completely frozen and the model maps onto the 
Jia-Jib-J2 model. Above T = O.lJi, the onset of orbital 
excitations causes a scrambling of the spin spectral weight, 
leading to sharply diminished spin-wave peaks, (iii) There is 
a sudden increase in the uniform magnetic susceptibility just 
below the phase transition. Above the transition, the uniform 
susceptibility continues to increase up to fairly high temper- 
atures, with a slope significantly higher than that in the Jia- 
•hb-'h or the J1-J2 model, (iv) The behavior of the specific 
heat and the order-parameter at the transition are very close to 
the corresponding Ising model on the same lattice, once the 
temperatures are scaled by the peak values. This suggests that 
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FIG. 2: (Color online) Plots for (a) the specific heat Cv, and (b) the 
uniform magnetic susceptibility Xm for the three models considered. 
The temperature T is expressed in terms of Ji (or Jia). Compared 
to the other two spin-only models, there is a sharp peak in Cv and a 
larger slope in Xm in the spin-orbital model. 



the transition is continuous and of second order, belonging to 
the universality class of the 2D Ising model. 

Fig. 2 shows the specific heat Cy and uniform magnetic 
susceptibility Xm for the spin-orbital model; for compari- 
son we also plot the same quantities for the spin-i Jia-Ju- 
J2 (with Jib = -O.lJia, J2 = 0.4Jia), and J1-J2 (with 
J2 — Ji) models. A main difference in Cy between the spin- 
orbital and the other two spin-only models is the sharp peak 
at T ~' 0.23 Ji, an indication of a phase transition. 

For an AF ordered ground state, Xm will grow as T in- 
creases from r=0, and then turn down at some characteristic 
temperature associated with short-range magnetic order. Xm 
at T =0 should have a finite value due to gapless excitations 
(Goldstone modes) intrinsic to each model in the thermody- 
namic limit. This is not captured in ED due to finite size ef- 
fects. Nonetheless, one expects the results to be qualitatively 
valid near the peak and quantitatively valid above it [18]. With 
our parameters, the energy to flip a spin in the AF ground state 
is approximately Ji + J/ + 2J2 ~ 2Ji. Hence the T = 
Xm should be comparable to that of an isotropic square-lattice 
Heisenberg model with the same Ji [19, 20]. One then ex- 
pects for both the spin-orbital and the Jia- Jib- J2 models an 
identical susceptibility below T = 0.1 Ji, with a magnitude of 
~ 0.05/ Ji. The sharp difference is the sudden increase in Xm 
between T = O.lJi and the phase transition ~ T = 0.2 Ji. 

A direct way to locate the orbital ordering transition tem- 
perature Tc is through the orbital Ising susceptibihty xi- 
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The sum on a is over the 2^^ Ising configurations, with 
the probability of the a*'* configuration. Nt = ni is the 
sum of the Ising variables on the lattice tied to the a*'* Ising 
configuration. According to the definition, xi is N/ 4 at T = 
0, and monotonically decreases to the configuration averaged 
value as T increases. The peak in dxi /dT is a measure of Tc 
which happens at ~ 0.23 Ji, as indicated by Fig. 3(a). 

We can define an orbital entropy S'orb = —jjTlaPot^'^Po', 
which approaches ln(2) at high temperature. On the other 
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FIG. 3: (Color online) (a) The orbital-Lsing Susceptibility xi and its 
derivative with respect to T in the inset, (b) The total system ^system 
and orbital S'orb entropy in the spin-orbital model. In the vicinity of 
the phase transition i? In 2 entropy is lost. 



hand, the total system entropy S'system incorporating both spin 
and orbital degrees of freedom (obtained by integrating Cv / T 
with respect to T) approaches ln(4) per site as T increases. 
Fig. 3(b) indicates that S'orb is completely exhausted soon after 
Tc, saturating much faster than 5system- We have checked that 
the behavior of the Ising variables in the spin-orbital model 
is quantitatively very close to the pure Ising model once the 
temperatures are scaled according to their corresponding Cv 
peak values. This suggests that the orbital phase transition is 
in the 2D Ising universality class where the finite temperature 
phase transition is continuous and of second order. A more 
definitive conclusion would require study via other numeri- 
cal techniques such as quantum Monte Carlo on larger sys- 
tems. This, however, may face minus sign problems because 
the spin-orbital model is frustrated. 

We next turn our focus to the spin dynamics of the spin- 
orbital model by studying the dynamic form factor S'"'^(q, w), 
which is the Fourier transform of the spin-spin correlation 
function {Sf{t) ■ S'j{t')). We calculate ^^^^(q,^^) via both 
ED and linear spin wave theory. Apart from the small energy 
gap in ED due to finite size effects, the spin wave dispersions 
obtained from both methods are compatible, see Fig. 4 (a). 
Neutron scattering on Fe-pnictide parent compounds indicates 
that the magnon energy is a maximum at (tt, tt). This behav- 
ior, absent in the J1-J2 model, is captured correctly by the 
Jia-Jib-Ji model, and hence the T = spin-orbtial model. 

At finite temperatures, orbital fluctuations start to play a 
role. The spectra in the spin-orbital model broaden much 
faster then the Jia-Jib-J2 model and show anomalous shifts 
to low frequencies. Fig. 4 (b) and (c) show the temperature 
evolution of spin wave at (tf, tt) obtained from ED. At temper- 
atures higher than 0.3 Ji, only incoherent spin waves sur- 
vive in the spin-orbital model. In contrast, in the J la- Jib- J 2 
model the coherent spin waves persist to a temperature higher 
than T ^ O.6J1. This feature can be seen also in the u- 
integrated form factors 5'"''(q) = J dujS°'^{q,uj). For the 
spin-orbital model, from T = to 0.4 Ji the dominant peak 
at (tt, 0) decreases by 25% in intensity, while for the Jia-Jib- 
J2 model it requires a temperature higher than T — 0.8 Ji to 
show a similar reduction, see Fig. 4 (d). Thus finite tempera- 



■■ Linear spin wave theory 

16-site Exact Diagonalization 



k (I 




((1,0)0 




2 4 

spin-wave energy a> 



^0.75 



CO 



0.5- 



(d) 


1 1 

— spin-orbital model 




— Jia-Jib-Ji model 







temperature T 



1.5 



HG. 4: (Color online) (a) T = (q, to) for the spin-orbital/ Jia- 
J16-J2 models, (b)-(c) Finite temperature S^'{(\,uj) at (7r,7r) ob- 
tained from ED for the spin-orbital, and the Jia- Jib- J2 models, re- 
spectively. The temperature goes from T = O.lJi (the blue curve) 
to T = l.OJi (the red curve), with a temperature increment between 
each curve ~ 0.08 Ji . A highly incoherent spin dynamics is observed 
in the spin-orbital model, (d) Finite temperature a;-integrated form 
factor S'^^(q) at (tt.O). 



ture neutron spectra can distinguish these models. 

Before we continue to discuss the relevance of this study 
to the Fe pnictides, a few comments are in order. The issue 
regarding the correlation strength in the Fe-pnictide materials 
is controversial and currently under debate. Recently, x-ray 
absorption data on several Fe-pnictide compounds revealed 
that the on-site Coulomb repulsion was smaller than the band- 
width, but also found a substantial Hund's coupling Jh = 0.8 
eV between the Fe 3d orbitals [21]. Moreover, there is no 
particular energy scale above which damped spin waves are 
found [4]; this absence of a Stoner decay strongly favors a 
picture based on localized moments. 

In many regards, the pnictides are schizophrenic, having 
aspects such as metallicity and strong covalency where cor- 
relations play a minor role [21, 22], and anti-ferromagnetism 
and local properties which derive directly from the strength 
of the Hund's coupling. Therefore a model based on local 
moments which takes aim at the magnetic properties of the 
pnictides and other transition metal oxide is completely in line 
with the findings in Ref. [21], and more recently with ob- 



4 



servations from optical conductivity measurements [23]. Our 
model focuses on a subset of localized orbitals in connection 
with magnetism, but neglects the fact that the 5 Fe d orbitals 
in Fe pnictides are not strongly crystal field split. 

Certain details of the model can be modified easily without 
changing the essential features. For example, there can be a 
direct coupling between the Ising variables reflecting lattice 
effects and quadrupolar couplings. In addition, the local en- 
vironment of As positions could modify local field screening 
and exchanges. These changes will alter the orbital gap and 
transition temperature, but not the overall picture. We will 
make our comparisons only in semi-quantitative terms. 

With this in mind, the spin-orbital model captures many 
features of the uniform susceptibility in the pnictides [24]. 
For example, in BaFe2As2 the susceptibility in emu/mole is 
0.6 X 10""^ at T = 0; it sharply increases near T = 150K 
to about 0.9 X 10^'^ and then continues to increase linearly 
to about 1.5 X 10"^ at T = 600A'. If part of the suscepti- 
bility is a weakly temperature dependent Pauli term, this im- 
plies an increase by a factor of about 3 between T = and 
T = 600K. Our finite-size calculations can be converted 
[25] to emu/mole by multiplying the susceptibility by a fac- 
tor {8* * 0.0938) / Ji, where Ji is in Kelvin and a factor of 
2 comes from the 2 Fe atoms per mole of the material. With 
Ji « lOOOX, this gives a susceptibility in cgs units of order 
lO^'^ at 500^, which comes down by about a factor of 3 — 4 
by T = including a sharp drop below the transition. 

In the 1111 pnictide family, two phase transitions at nearby 
temperatures have been reported, a structural transition at 
higher T and a magnetic transition at lower T. In contrast, 
only one simultaneous structural and magnetic transition is 
found in the 122 family. This is explained naturally in terms 
of 3D couplings. Orbital ordering driven by magnetism re- 
quires the prior development of short-range spin order. A 
3D system with strong inter-planar coupling would therefore 
lead to simultaneous spin and orbital order In contrast, for 
a weakly inter-planar coupled 2D system orbitals order when 
short-range spin order develops, but spins only order on the 
scale of inter-planar couplings. This therefore leads to two 
separate transitions. These observations are indeed consistent 
with the two different pnictide families. This aspect also has 
been suggested for J1-J2 Heisenberg models [6]. However, 
one important difference is that in the spin-orbital model ~ 
R In 2 entropy is lost in the vicinity of the transition; it is likely 
significantly smaller in the J1-J2 models. It is noted that in 
the Fei+ySe^^Tei-j; systems a comparable amount of entropy 
change is found near the AF transition [26] . 

In summary, we have studied a model system that captures 
the physics of coupled spin and orbital degrees of freedom. 
Such a system apparently undergoes a continuous, second- 
order phase transition to an orbitally ordered state at a tem- 
perature set by short-range magnetic order The onset of or- 
bital excitations and fluctuations cause a highly incoherent 
spin dynamics, leading to a sharp increase in uniform mag- 
netic susceptibility. The susceptibility continues to increase 
up to fairly high temperature above the phase transition, with 



a large slope comparable to those observed in the pnictides. 
Our calculations of dynamic structure factors at finite tem- 
peratures serve as clear predictions of the spin-orbital model, 
that can be tested by further experiments. In addition to the 
pnictides, the model should be generally applicable to other 
systems with orbital degeneracy, with the strengths/signs of 
the exchange constants dependent on the microscopic details. 
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